In [1]:
import mesa as ms
# quantidade de zonas no modelo 500
folder2 = '/home/glauffer/Dropbox/UFRGS/kepler/mesa/alpha2/'
a = ms.mesa_profile(folder2, 500, num_type='nearest_model')
58 in profiles.index file ...
Found and load nearest profile for cycle 500
reading /home/glauffer/Dropbox/UFRGS/kepler/mesa/alpha2//profile16.data ...
reading ...100%
In [2]:
a.cols
Out[2]:
{'abar': 29,
'am_log_D_DSI': 117,
'am_log_D_ES': 120,
'am_log_D_GSF': 121,
'am_log_D_SH': 118,
'am_log_D_SSI': 119,
'am_log_D_ST': 122,
'am_log_nu_ST': 123,
'am_log_sig': 116,
'burn_ar': 86,
'burn_c': 78,
'burn_ca': 87,
'burn_cr': 89,
'burn_fe': 90,
'burn_mg': 83,
'burn_n': 79,
'burn_na': 82,
'burn_ne': 81,
'burn_o': 80,
'burn_s': 85,
'burn_si': 84,
'burn_ti': 88,
'c12': 62,
'c12_c12': 91,
'c12_o16': 92,
'c13': 63,
'cno': 76,
'conv_vel_div_csound': 44,
'csound': 11,
'dlnT_dt': 55,
'dlnd_dt': 54,
'dq_ratio': 18,
'dynamo_log_B_phi': 125,
'dynamo_log_B_r': 124,
'entropy': 9,
'eps_grav': 34,
'eps_nuc': 32,
'eta': 15,
'f19': 70,
'fp_rot': 106,
'free_e': 28,
'ft_rot': 107,
'gradL': 51,
'gradT': 47,
'gradT_sub_grada': 37,
'grad_density': 50,
'grad_temperature': 49,
'grada': 27,
'gradr': 48,
'gradr_sub_grada': 126,
'h1': 59,
'he3': 60,
'he4': 61,
'i_rot': 104,
'j_rot': 105,
'ledoux_stable': 53,
'logL': 7,
'logP': 4,
'logPgas': 26,
'logR': 5,
'logRho': 3,
'logT': 2,
'log_D_conv': 38,
'log_D_mix': 42,
'log_D_ovr': 40,
'log_D_semi': 39,
'log_D_thrm': 41,
'log_J_div_M53': 100,
'log_J_inside': 101,
'log_abs_shear': 102,
'log_am_nu': 108,
'log_conv_vel': 43,
'log_j_rot': 99,
'log_mlt_D_mix': 45,
'log_omega': 98,
'log_opacity': 31,
'log_r_equatorial': 112,
'log_r_polar': 110,
'logdq': 17,
'logtau': 23,
'logxq': 58,
'luminosity': 6,
'mass': 56,
'mg24': 73,
'mixing_type': 10,
'mlt_mixing_length': 35,
'mlt_mixing_type': 36,
'mmid': 57,
'mu': 16,
'n13': 64,
'n14': 65,
'n15': 66,
'ne20': 71,
'ne22': 72,
'non_nuc_neu': 33,
'o16': 67,
'o16_o16': 93,
'o17': 68,
'o18': 69,
'omega': 97,
'omega_crit': 114,
'omega_div_omega_crit': 115,
'other': 96,
'pgas_div_ptotal': 25,
'photo': 95,
'pnhe4': 94,
'pp': 75,
'pressure': 24,
'pressure_scale_height': 46,
'q': 19,
'r_e_div_r_p': 113,
'r_equatorial': 111,
'r_polar': 109,
'radius': 20,
'sch_stable': 52,
'si28': 74,
'tau': 22,
'temperature': 21,
'total_energy': 13,
'tri_alfa': 77,
'v_div_csound': 12,
'v_div_r': 14,
'v_rot': 103,
'velocity': 8,
'ye': 30,
'zone': 1}
In [3]:
z = a.get('zone')
In [4]:
z
Out[4]:
array([ 1., 2., 3., 4., 5., 6., 7., 8., 9.,
10., 11., 12., 13., 14., 15., 16., 17., 18.,
19., 20., 21., 22., 23., 24., 25., 26., 27.,
28., 29., 30., 31., 32., 33., 34., 35., 36.,
37., 38., 39., 40., 41., 42., 43., 44., 45.,
46., 47., 48., 49., 50., 51., 52., 53., 54.,
55., 56., 57., 58., 59., 60., 61., 62., 63.,
64., 65., 66., 67., 68., 69., 70., 71., 72.,
73., 74., 75., 76., 77., 78., 79., 80., 81.,
82., 83., 84., 85., 86., 87., 88., 89., 90.,
91., 92., 93., 94., 95., 96., 97., 98., 99.,
100., 101., 102., 103., 104., 105., 106., 107., 108.,
109., 110., 111., 112., 113., 114., 115., 116., 117.,
118., 119., 120., 121., 122., 123., 124., 125., 126.,
127., 128., 129., 130., 131., 132., 133., 134., 135.,
136., 137., 138., 139., 140., 141., 142., 143., 144.,
145., 146., 147., 148., 149., 150., 151., 152., 153.,
154., 155., 156., 157., 158., 159., 160., 161., 162.,
163., 164., 165., 166., 167., 168., 169., 170., 171.,
172., 173., 174., 175., 176., 177., 178., 179., 180.,
181., 182., 183., 184., 185., 186., 187., 188., 189.,
190., 191., 192., 193., 194., 195., 196., 197., 198.,
199., 200., 201., 202., 203., 204., 205., 206., 207.,
208., 209., 210., 211., 212., 213., 214., 215., 216.,
217., 218., 219., 220., 221., 222., 223., 224., 225.,
226., 227., 228., 229., 230., 231., 232., 233., 234.,
235., 236., 237., 238., 239., 240., 241., 242., 243.,
244., 245., 246., 247., 248., 249., 250., 251., 252.,
253., 254., 255., 256., 257., 258., 259., 260., 261.,
262., 263., 264., 265., 266., 267., 268., 269., 270.,
271., 272., 273., 274., 275., 276., 277., 278., 279.,
280., 281., 282., 283., 284., 285., 286., 287., 288.,
289., 290., 291., 292., 293., 294., 295., 296., 297.,
298., 299., 300., 301., 302., 303., 304., 305., 306.,
307., 308., 309., 310., 311., 312., 313., 314., 315.,
316., 317., 318., 319., 320., 321., 322., 323., 324.,
325., 326., 327., 328., 329., 330., 331., 332., 333.,
334., 335., 336., 337., 338., 339., 340., 341., 342.,
343., 344., 345., 346., 347., 348., 349., 350., 351.,
352., 353., 354., 355., 356., 357., 358., 359., 360.,
361., 362., 363., 364., 365., 366., 367., 368., 369.,
370., 371., 372., 373., 374., 375., 376., 377., 378.,
379., 380., 381., 382., 383., 384., 385., 386., 387.,
388., 389., 390., 391., 392., 393., 394., 395., 396.,
397., 398., 399., 400., 401., 402., 403., 404., 405.,
406., 407., 408., 409., 410., 411., 412., 413., 414.,
415., 416., 417., 418., 419., 420., 421., 422., 423.,
424., 425., 426., 427., 428., 429., 430., 431., 432.,
433., 434., 435., 436., 437., 438., 439., 440., 441.,
442., 443., 444., 445., 446., 447., 448., 449., 450.,
451., 452., 453., 454., 455., 456., 457., 458., 459.,
460., 461., 462., 463., 464., 465., 466., 467., 468.,
469., 470., 471., 472., 473., 474., 475., 476., 477.,
478., 479., 480., 481., 482., 483., 484., 485., 486.,
487., 488., 489., 490., 491., 492., 493., 494., 495.,
496., 497., 498., 499., 500., 501., 502., 503., 504.,
505., 506., 507., 508., 509., 510., 511., 512., 513.,
514., 515., 516., 517., 518., 519., 520., 521., 522.,
523., 524., 525., 526., 527., 528., 529., 530., 531.,
532., 533., 534., 535., 536., 537., 538., 539., 540.,
541., 542., 543., 544., 545., 546., 547., 548., 549.,
550., 551., 552., 553., 554., 555., 556., 557., 558.,
559., 560., 561., 562., 563., 564., 565., 566., 567.,
568., 569., 570., 571., 572., 573., 574., 575., 576.,
577., 578., 579., 580., 581., 582., 583., 584., 585.,
586., 587., 588., 589., 590., 591., 592., 593., 594.,
595., 596., 597., 598., 599., 600., 601., 602., 603.,
604., 605., 606., 607., 608., 609., 610., 611., 612.,
613., 614., 615., 616., 617., 618., 619., 620., 621.,
622., 623., 624., 625., 626., 627., 628., 629., 630.,
631., 632., 633., 634., 635., 636., 637., 638., 639.,
640., 641., 642., 643., 644., 645., 646., 647., 648.,
649., 650., 651., 652., 653., 654., 655., 656., 657.,
658., 659., 660., 661., 662., 663., 664., 665., 666.,
667., 668., 669., 670., 671., 672., 673., 674., 675.,
676., 677., 678., 679., 680., 681., 682., 683., 684.,
685., 686., 687., 688., 689., 690., 691., 692., 693.,
694., 695., 696., 697., 698., 699., 700., 701., 702.,
703., 704., 705., 706., 707., 708., 709., 710., 711.,
712., 713., 714., 715., 716., 717., 718., 719., 720.,
721., 722., 723., 724., 725., 726., 727., 728., 729.,
730., 731., 732., 733., 734., 735., 736., 737., 738.,
739., 740., 741., 742., 743., 744., 745., 746., 747.,
748., 749., 750., 751., 752., 753., 754., 755., 756.,
757., 758., 759., 760., 761., 762., 763., 764., 765.,
766., 767., 768., 769., 770., 771., 772., 773., 774.,
775., 776., 777., 778., 779., 780., 781., 782., 783.,
784., 785., 786., 787., 788., 789., 790., 791., 792.,
793., 794., 795., 796., 797., 798., 799., 800., 801.,
802., 803., 804., 805., 806., 807., 808., 809., 810.,
811., 812., 813., 814., 815., 816., 817., 818., 819.,
820., 821., 822., 823., 824., 825., 826., 827., 828.,
829., 830., 831., 832., 833., 834., 835., 836., 837.,
838., 839., 840., 841., 842., 843., 844., 845., 846.,
847., 848., 849., 850., 851., 852., 853., 854., 855.,
856., 857., 858., 859., 860., 861., 862., 863., 864.,
865., 866., 867., 868., 869., 870., 871., 872., 873.,
874., 875., 876., 877., 878., 879., 880., 881., 882.,
883., 884., 885., 886., 887., 888., 889., 890., 891.,
892., 893., 894., 895., 896., 897., 898., 899., 900.,
901., 902., 903., 904., 905., 906., 907., 908., 909.,
910., 911., 912., 913., 914., 915., 916., 917., 918.,
919., 920., 921., 922., 923., 924., 925., 926., 927.,
928., 929., 930., 931., 932., 933., 934., 935., 936.,
937., 938., 939., 940., 941., 942., 943., 944., 945.,
946., 947., 948., 949., 950., 951., 952., 953., 954.,
955., 956., 957., 958., 959., 960., 961., 962., 963.,
964., 965., 966., 967., 968., 969., 970., 971., 972.,
973., 974., 975., 976., 977., 978., 979., 980., 981.,
982., 983., 984., 985., 986., 987., 988., 989., 990.,
991., 992., 993., 994.])
In [6]:
# quantidade de zonas no modelo 2090
a2090 = ms.mesa_profile(folder2, 2090, num_type='nearest_model')
z2090 = a2090.get('zone')
z2090
58 in profiles.index file ...
Found and load nearest profile for cycle 2100
reading /home/glauffer/Dropbox/UFRGS/kepler/mesa/alpha2//profile48.data ...
reading ...100%
Out[6]:
array([ 1., 2., 3., 4., 5., 6., 7., 8., 9.,
10., 11., 12., 13., 14., 15., 16., 17., 18.,
19., 20., 21., 22., 23., 24., 25., 26., 27.,
28., 29., 30., 31., 32., 33., 34., 35., 36.,
37., 38., 39., 40., 41., 42., 43., 44., 45.,
46., 47., 48., 49., 50., 51., 52., 53., 54.,
55., 56., 57., 58., 59., 60., 61., 62., 63.,
64., 65., 66., 67., 68., 69., 70., 71., 72.,
73., 74., 75., 76., 77., 78., 79., 80., 81.,
82., 83., 84., 85., 86., 87., 88., 89., 90.,
91., 92., 93., 94., 95., 96., 97., 98., 99.,
100., 101., 102., 103., 104., 105., 106., 107., 108.,
109., 110., 111., 112., 113., 114., 115., 116., 117.,
118., 119., 120., 121., 122., 123., 124., 125., 126.,
127., 128., 129., 130., 131., 132., 133., 134., 135.,
136., 137., 138., 139., 140., 141., 142., 143., 144.,
145., 146., 147., 148., 149., 150., 151., 152., 153.,
154., 155., 156., 157., 158., 159., 160., 161., 162.,
163., 164., 165., 166., 167., 168., 169., 170., 171.,
172., 173., 174., 175., 176., 177., 178., 179., 180.,
181., 182., 183., 184., 185., 186., 187., 188., 189.,
190., 191., 192., 193., 194., 195., 196., 197., 198.,
199., 200., 201., 202., 203., 204., 205., 206., 207.,
208., 209., 210., 211., 212., 213., 214., 215., 216.,
217., 218., 219., 220., 221., 222., 223., 224., 225.,
226., 227., 228., 229., 230., 231., 232., 233., 234.,
235., 236., 237., 238., 239., 240., 241., 242., 243.,
244., 245., 246., 247., 248., 249., 250., 251., 252.,
253., 254., 255., 256., 257., 258., 259., 260., 261.,
262., 263., 264., 265., 266., 267., 268., 269., 270.,
271., 272., 273., 274., 275., 276., 277., 278., 279.,
280., 281., 282., 283., 284., 285., 286., 287., 288.,
289., 290., 291., 292., 293., 294., 295., 296., 297.,
298., 299., 300., 301., 302., 303., 304., 305., 306.,
307., 308., 309., 310., 311., 312., 313., 314., 315.,
316., 317., 318., 319., 320., 321., 322., 323., 324.,
325., 326., 327., 328., 329., 330., 331., 332., 333.,
334., 335., 336., 337., 338., 339., 340., 341., 342.,
343., 344., 345., 346., 347., 348., 349., 350., 351.,
352., 353., 354., 355., 356., 357., 358., 359., 360.,
361., 362., 363., 364., 365., 366., 367., 368., 369.,
370., 371., 372., 373., 374., 375., 376., 377., 378.,
379., 380., 381., 382., 383., 384., 385., 386., 387.,
388., 389., 390., 391., 392., 393., 394., 395., 396.,
397., 398., 399., 400., 401., 402., 403., 404., 405.,
406., 407., 408., 409., 410., 411., 412., 413., 414.,
415., 416., 417., 418., 419., 420., 421., 422., 423.,
424., 425., 426., 427., 428., 429., 430., 431., 432.,
433., 434., 435., 436., 437., 438., 439., 440., 441.,
442., 443., 444., 445., 446., 447., 448., 449., 450.,
451., 452., 453., 454., 455., 456., 457., 458., 459.,
460., 461., 462., 463., 464., 465., 466., 467., 468.,
469., 470., 471., 472., 473., 474., 475., 476., 477.,
478., 479., 480., 481., 482., 483., 484., 485., 486.,
487., 488., 489., 490., 491., 492., 493., 494., 495.,
496., 497., 498., 499., 500., 501., 502., 503., 504.,
505., 506., 507., 508., 509., 510., 511., 512., 513.,
514., 515., 516., 517., 518., 519., 520., 521., 522.,
523., 524., 525., 526., 527., 528., 529., 530., 531.,
532., 533., 534., 535., 536., 537., 538., 539., 540.,
541., 542., 543., 544., 545., 546., 547., 548., 549.,
550., 551., 552., 553., 554., 555., 556., 557., 558.,
559., 560., 561., 562., 563., 564., 565., 566., 567.,
568., 569., 570., 571., 572., 573., 574., 575., 576.,
577., 578., 579., 580., 581., 582., 583., 584., 585.,
586., 587., 588., 589., 590., 591., 592., 593., 594.,
595., 596., 597., 598., 599., 600., 601., 602., 603.,
604., 605., 606., 607., 608., 609., 610., 611., 612.,
613., 614., 615., 616., 617., 618., 619., 620., 621.,
622., 623., 624., 625., 626., 627., 628., 629., 630.,
631., 632., 633., 634., 635., 636., 637., 638., 639.,
640., 641., 642., 643., 644., 645., 646., 647., 648.,
649., 650., 651., 652., 653., 654., 655., 656., 657.,
658., 659., 660., 661., 662., 663., 664., 665., 666.,
667., 668., 669., 670., 671., 672., 673., 674., 675.,
676., 677., 678., 679., 680., 681., 682., 683., 684.,
685., 686., 687., 688., 689., 690., 691., 692., 693.,
694., 695., 696., 697., 698., 699., 700., 701., 702.,
703., 704., 705., 706., 707., 708., 709., 710., 711.,
712., 713., 714., 715., 716., 717., 718., 719., 720.,
721., 722., 723., 724., 725., 726., 727., 728., 729.,
730., 731., 732., 733., 734., 735., 736., 737., 738.,
739., 740., 741., 742., 743., 744., 745., 746., 747.,
748., 749., 750., 751., 752., 753., 754., 755., 756.,
757., 758., 759., 760., 761., 762., 763., 764., 765.,
766., 767., 768., 769., 770., 771., 772., 773., 774.,
775., 776., 777., 778., 779., 780., 781., 782., 783.,
784., 785., 786., 787., 788., 789., 790., 791., 792.,
793., 794., 795., 796., 797., 798., 799., 800., 801.,
802., 803., 804., 805., 806., 807., 808., 809., 810.,
811., 812., 813., 814., 815., 816., 817., 818., 819.,
820., 821., 822., 823., 824., 825., 826., 827., 828.,
829., 830., 831., 832., 833., 834., 835., 836., 837.,
838., 839., 840., 841., 842., 843., 844., 845., 846.,
847., 848., 849., 850., 851., 852., 853., 854., 855.,
856., 857., 858., 859., 860., 861., 862., 863., 864.,
865., 866., 867., 868., 869., 870., 871., 872., 873.,
874., 875., 876., 877., 878., 879., 880., 881., 882.,
883., 884., 885., 886., 887., 888., 889., 890., 891.,
892., 893., 894., 895., 896., 897., 898., 899., 900.,
901., 902., 903., 904., 905., 906., 907., 908., 909.,
910., 911., 912., 913., 914., 915., 916., 917., 918.,
919., 920., 921., 922., 923., 924., 925., 926., 927.,
928., 929., 930., 931., 932., 933., 934., 935., 936.,
937., 938., 939., 940., 941., 942., 943., 944., 945.,
946., 947., 948., 949., 950., 951., 952., 953., 954.,
955., 956., 957., 958., 959., 960., 961., 962., 963.,
964., 965., 966., 967., 968., 969.])
In [8]:
t2090 = a2090.get('logT')
In [9]:
t2090
Out[9]:
array([ 3.78719491, 3.78719492, 3.78719493, 3.78719495, 3.78719496,
3.78719499, 3.78719503, 3.78719509, 3.78719516, 3.78719525,
3.78719536, 3.78719554, 3.7871958 , 3.78719611, 3.78719649,
3.78719704, 3.78719775, 3.78719894, 3.78720058, 3.78720294,
3.787206 , 3.78721017, 3.78721544, 3.78722275, 3.78723212,
3.78724485, 3.78726095, 3.78728216, 3.78730851, 3.78734232,
3.78738365, 3.78743523, 3.78749714, 3.78757283, 3.78766243,
3.78777307, 3.78790501, 3.7880742 , 3.78828133, 3.78855846,
3.78890771, 3.78939223, 3.79001915, 3.79089798, 3.79205266,
3.79363935, 3.79573115, 3.7985719 , 3.80218789, 3.80680435,
3.81243658, 3.82021253, 3.83063209, 3.84737008, 3.87571331,
3.91077412, 3.94239392, 3.9590003 , 3.97420115, 3.98988484,
4.00849633, 4.0335859 , 4.05981097, 4.07882079, 4.09312834,
4.10457603, 4.11414465, 4.1223866 , 4.12963993, 4.13613103,
4.14202205, 4.14743291, 4.15158667, 4.15472422, 4.15773667,
4.16063751, 4.16343825, 4.16614879, 4.16877767, 4.17133231,
4.1738192 , 4.17624405, 4.17979131, 4.18432609, 4.18868756,
4.19445698, 4.20141835, 4.20805629, 4.21441495, 4.22052493,
4.22640821, 4.23208141, 4.24026815, 4.25048489, 4.26004018,
4.26899897, 4.27866584, 4.28996337, 4.29859173, 4.30410354,
4.31204103, 4.3207097 , 4.32777506, 4.33340573, 4.33994392,
4.3482184 , 4.35562725, 4.36327802, 4.37287358, 4.38064963,
4.38816457, 4.39651374, 4.40339965, 4.40894661, 4.41543073,
4.42373638, 4.43120717, 4.43743934, 4.44440392, 4.4537991 ,
4.46406688, 4.47389712, 4.48332449, 4.49013269, 4.49677337,
4.50590158, 4.51635712, 4.52796339, 4.53813279, 4.54517266,
4.55069148, 4.55746565, 4.56675056, 4.5771423 , 4.58686957,
4.59448904, 4.60053824, 4.60733503, 4.61484357, 4.62296388,
4.63307797, 4.64509223, 4.65650487, 4.66496642, 4.67324551,
4.68453845, 4.69674507, 4.70529312, 4.71070653, 4.71607252,
4.72138055, 4.72792904, 4.73567148, 4.74203314, 4.74706753,
4.75453318, 4.76435872, 4.77266723, 4.77947108, 4.78681711,
4.79660864, 4.80618392, 4.81186067, 4.81561493, 4.82120165,
4.82858709, 4.83965971, 4.8507882 , 4.85821651, 4.86466841,
4.87015833, 4.87803805, 4.88573523, 4.8908242 , 4.89818774,
4.90779258, 4.91617667, 4.92335857, 4.93166687, 4.93963371,
4.94607878, 4.95335697, 4.96146803, 4.96900416, 4.97596983,
4.98383558, 4.99260894, 5.00063302, 5.00788574, 5.01530633,
5.02289132, 5.03074098, 5.03684122, 5.04088531, 5.04584338,
5.05171253, 5.06019856, 5.06927575, 5.07626055, 5.08344753,
5.09083824, 5.09842874, 5.10607307, 5.11322981, 5.12005093,
5.12807656, 5.13731549, 5.14499742, 5.15111339, 5.16023824,
5.17007832, 5.17752113, 5.18495386, 5.19237851, 5.19984273,
5.20734899, 5.21485255, 5.22151748, 5.22912721, 5.23661172,
5.24210783, 5.24721696, 5.25193416, 5.2580185 , 5.26549551,
5.27230146, 5.27842274, 5.28336687, 5.28711467, 5.2919323 ,
5.29783809, 5.30290742, 5.30712513, 5.31274788, 5.3198061 ,
5.32698432, 5.3324727 , 5.33855058, 5.3473101 , 5.35386636,
5.35835414, 5.36325801, 5.3706116 , 5.37804673, 5.38300289,
5.38758666, 5.39179018, 5.39739949, 5.40445199, 5.41051907,
5.41557169, 5.42001007, 5.42382084, 5.42889805, 5.4363717 ,
5.44477437, 5.4529798 , 5.46002696, 5.46576742, 5.47230524,
5.47801775, 5.48207639, 5.4870916 , 5.49303587, 5.49826301,
5.50279638, 5.50761447, 5.5127093 , 5.51754976, 5.52214075,
5.52690041, 5.53182515, 5.53630306, 5.5403428 , 5.54475168,
5.54951986, 5.5542433 , 5.55832701, 5.56178517, 5.56578608,
5.56989192, 5.57499026, 5.57966888, 5.58287751, 5.58641233,
5.59165562, 5.5971668 , 5.60124834, 5.60528377, 5.60927294,
5.6131152 , 5.61696348, 5.62052022, 5.62403856, 5.62790073,
5.63133477, 5.63510693, 5.63882876, 5.64214062, 5.64541385,
5.64900996, 5.6543213 , 5.65952887, 5.6636297 , 5.6670003 ,
5.6709966 , 5.67590562, 5.68017003, 5.68452737, 5.68900695,
5.69398616, 5.69924234, 5.70384563, 5.7087963 , 5.71418588,
5.71940938, 5.72421913, 5.72863652, 5.73267692, 5.736967 ,
5.74088976, 5.74476315, 5.74863841, 5.75292643, 5.75762839,
5.76180208, 5.76547308, 5.76910011, 5.77445542, 5.78005269,
5.78520591, 5.79060194, 5.79622856, 5.80061699, 5.80398923,
5.80992702, 5.81620727, 5.82192674, 5.8275335 , 5.83194795,
5.83629896, 5.84058833, 5.8453434 , 5.85105499, 5.85564442,
5.8596711 , 5.86413391, 5.86853355, 5.8723883 , 5.87572372,
5.87973997, 5.88536325, 5.89111804, 5.89452003, 5.89676847,
5.90011216, 5.9034162 , 5.90668837, 5.91099519, 5.9162977 ,
5.92253999, 5.9286606 , 5.93367086, 5.9376252 , 5.94153155,
5.94539114, 5.95125939, 5.95943989, 5.96574559, 5.96988551,
5.97600224, 5.98215656, 5.98679795, 5.99222476, 5.99865706,
6.00558776, 6.01237819, 6.01903397, 6.0255606 , 6.0319635 ,
6.03824791, 6.04480158, 6.05160681, 6.06157986, 6.07126596,
6.07758847, 6.08380492, 6.08992077, 6.09792915, 6.1077154 ,
6.117275 , 6.1243064 , 6.12893138, 6.13577724, 6.14474068,
6.15352528, 6.16213836, 6.16975001, 6.17640877, 6.18297082,
6.18944028, 6.19897602, 6.21140651, 6.220537 , 6.22653162,
6.23538809, 6.24696682, 6.25548682, 6.26108531, 6.26662065,
6.2720932 , 6.28018189, 6.29074822, 6.29851063, 6.30360784,
6.31113607, 6.32095077, 6.3305228 , 6.33986169, 6.34671878,
6.35122532, 6.35789011, 6.36660118, 6.37512646, 6.38347515,
6.39165516, 6.39967329, 6.40753554, 6.41715293, 6.42561481,
6.43116126, 6.43844044, 6.44735732, 6.4560842 , 6.46463058,
6.47134442, 6.47631138, 6.4836552 , 6.49246762, 6.50033366,
6.50806986, 6.51568408, 6.5239261 , 6.5327671 , 6.5407389 ,
6.54610435, 6.54964882, 6.55491713, 6.56185692, 6.56937879,
6.57745599, 6.58474611, 6.59128048, 6.59772933, 6.60409528,
6.61224859, 6.62210762, 6.63058208, 6.63773409, 6.64711237,
6.65575939, 6.66144119, 6.66984176, 6.68083222, 6.68945559,
6.69582248, 6.70210516, 6.70830447, 6.71442134, 6.72045692,
6.72935919, 6.73808847, 6.74381333, 6.75226196, 6.76055057,
6.76598973, 6.7740212 , 6.7843529 , 6.79429249, 6.80124345,
6.80535754, 6.81144887, 6.81941722, 6.82721711, 6.83485214,
6.84232571, 6.84964132, 6.85680197, 6.86438898, 6.87236023,
6.88067931, 6.88929967, 6.89662896, 6.90276352, 6.9087605 ,
6.91510741, 6.9217668 , 6.9310913 , 6.94097536, 6.94795586,
6.95391528, 6.96258967, 6.97087961, 6.97891555, 6.98775823,
6.99465582, 7.00129306, 7.01013979, 7.01865516, 7.02487347,
7.03183662, 7.03824452, 7.04460233, 7.05232023, 7.05959162,
7.06964763, 7.07947473, 7.08668327, 7.09543542, 7.10284034,
7.10825546, 7.11486983, 7.12261384, 7.13035983, 7.13643637,
7.14095483, 7.14761073, 7.15738663, 7.16801183, 7.17642127,
7.18378169, 7.19298257, 7.20156806, 7.20864627, 7.21555645,
7.22443154, 7.23467084, 7.24219615, 7.24784547, 7.255633 ,
7.26542871, 7.27303778, 7.27836645, 7.28594347, 7.29824791,
7.30922296, 7.31511397, 7.32250458, 7.33342139, 7.34138407,
7.34585954, 7.35092084, 7.35654208, 7.36372663, 7.37254393,
7.38022687, 7.38668759, 7.39316309, 7.40051113, 7.40990219,
7.41929807, 7.42688938, 7.43508206, 7.44313496, 7.45031635,
7.4578655 , 7.46746592, 7.47756634, 7.48370706, 7.4874483 ,
7.492513 , 7.49886981, 7.50850097, 7.51776678, 7.52316824,
7.53172334, 7.54205446, 7.55073389, 7.55851694, 7.56582171,
7.57566918, 7.58528857, 7.59247373, 7.60003667, 7.60688611,
7.61259656, 7.6175382 , 7.6217182 , 7.62755095, 7.6350236 ,
7.6420222 , 7.65049319, 7.65885893, 7.6651919 , 7.67153369,
7.68223233, 7.69354102, 7.70096636, 7.71119433, 7.72166419,
7.72946345, 7.73853651, 7.74588051, 7.75014497, 7.75600104,
7.76346947, 7.77086631, 7.77819326, 7.78611127, 7.79463718,
7.80208552, 7.80843173, 7.81600824, 7.8240302 , 7.831236 ,
7.84121128, 7.85392944, 7.8634568 , 7.86981198, 7.88087084,
7.8924036 , 7.89932809, 7.9092709 , 7.92143414, 7.92991895,
7.93559293, 7.94475106, 7.95626866, 7.96701489, 7.97561174,
7.9826188 , 7.99063272, 7.99836523, 8.00456623, 8.0088216 ,
8.01178358, 8.01376214, 8.0157469 , 8.01773797, 8.0197354 ,
8.02173929, 8.0237497 , 8.02576673, 8.02779046, 8.02982096,
8.03185833, 8.03492121, 8.03789725, 8.03973018, 8.04156886,
8.04341335, 8.04618594, 8.04898243, 8.05085083, 8.05272537,
8.05460614, 8.05649319, 8.0583866 , 8.06123326, 8.06410584,
8.06554633, 8.06650781, 8.06800561, 8.0700475 , 8.07209709,
8.07415447, 8.07621973, 8.07829298, 8.07985492, 8.08089762,
8.08246381, 8.08456162, 8.08666781, 8.0887825 , 8.09090579,
8.09303779, 8.0951786 , 8.09732835, 8.09948713, 8.10168329,
8.10391757, 8.10616173, 8.1084159 , 8.1106802 , 8.11295477,
8.11466961, 8.1158147 , 8.11753516, 8.11984126, 8.12215811,
8.12448585, 8.12682458, 8.12917442, 8.13153549, 8.1339079 ,
8.13629174, 8.1387879 , 8.14139911, 8.14336936, 8.14468532,
8.14666295, 8.14931564, 8.15198222, 8.15399415, 8.15533789,
8.15735719, 8.16026896, 8.16262241, 8.16419466, 8.16655798,
8.16893982, 8.170531 , 8.17212701, 8.17372799, 8.17613483,
8.17856136, 8.18018258, 8.18261968, 8.18507627, 8.18671733,
8.18836318, 8.19001382, 8.19249455, 8.1949937 , 8.19666251,
8.19916954, 8.20169285, 8.20337683, 8.205905 , 8.20844581,
8.21013997, 8.21175651, 8.21329443, 8.21559966, 8.21867434,
8.22097809, 8.2225103 , 8.22480231, 8.22708702, 8.22860414,
8.23086992, 8.23312073, 8.23461266, 8.23668532, 8.23932082,
8.24192418, 8.2444921 , 8.24702134, 8.24950965, 8.25195603,
8.25435958, 8.25671431, 8.25901951, 8.26240484, 8.26566247,
8.26778669, 8.26986682, 8.27190402, 8.27389947, 8.2758541 ,
8.27729094, 8.27823851, 8.27964456, 8.28148219, 8.28283316,
8.28372416, 8.28504638, 8.28677468, 8.28804336, 8.28887736,
8.29011525, 8.29173361, 8.29331959, 8.29487371, 8.29639647,
8.29788835, 8.29934982, 8.3007813 , 8.30218324, 8.30355605,
8.30490014, 8.30621595, 8.30750388, 8.30876436, 8.3101328 ,
8.31159725, 8.31273124, 8.31355933, 8.31478151, 8.31636141,
8.31789257, 8.31937688, 8.32081633, 8.32221309, 8.32325391,
8.32396743, 8.32502141, 8.32638735, 8.32738477, 8.32804098,
8.32868917, 8.3293298 , 8.32996331, 8.33059019, 8.33105579,
8.33136468, 8.33182584, 8.33243557, 8.33304051, 8.33357566,
8.33404284, 8.33474015, 8.33543218, 8.33589182, 8.33635009,
8.33680715, 8.33726319, 8.33771839, 8.33817291, 8.33862691,
8.33908054, 8.33953396, 8.34021379, 8.34089431, 8.34134815,
8.34202933, 8.34271246, 8.34316842, 8.34362511, 8.34408261,
8.34476996, 8.34569042, 8.34638401, 8.34684734, 8.34748003,
8.34828455, 8.34902814, 8.34970973, 8.35056454, 8.35133828,
8.35185527, 8.3524599 , 8.35315326, 8.35384888, 8.35437235,
8.35472179, 8.35524666, 8.35594872, 8.35665281, 8.35735887,
8.35788996, 8.35824442, 8.35877669, 8.3598268 , 8.36087417,
8.36157368, 8.36239141, 8.36332863, 8.36426826, 8.36497488,
8.36544639, 8.36615425, 8.36725785, 8.36820684, 8.36884006,
8.36978359, 8.37166697, 8.37473933, 8.37759139, 8.37981087,
8.38298341, 8.38656088, 8.38933979, 8.39187193, 8.39416327,
8.39731269, 8.40098673, 8.40365493, 8.4064646 , 8.40918854,
8.41171228, 8.41509272, 8.41819752, 8.42075984, 8.42416452,
8.4269814 , 8.42977954, 8.43283413, 8.43575787, 8.43965928,
8.44366426, 8.44692593, 8.44996215, 8.45260319, 8.45552121,
8.45964085, 8.46329958, 8.4671525 , 8.47136216, 8.47446316,
8.47731144, 8.47990979, 8.48266523, 8.48654149, 8.49005795,
8.49275728, 8.49564274, 8.49933896, 8.50327706, 8.50670804,
8.5102793 , 8.51345518, 8.51658413, 8.51951518, 8.52140966,
8.52415193, 8.5282429 , 8.53223711, 8.53482789, 8.53739072,
8.54093134, 8.54404029, 8.5461984 , 8.54874319, 8.55150881,
8.55419245, 8.55655847, 8.55917065, 8.56216913, 8.56399057,
8.56610841, 8.56914105, 8.57044859, 8.57337184, 8.57614926,
8.57846313, 8.58116662, 8.58365063, 8.58551063, 8.58724599,
8.5893088 , 8.59121019, 8.59319254, 8.59497477, 8.59582172,
8.5963644 , 8.59715853, 8.59775269, 8.59812504, 8.59870478,
8.5994718 , 8.6001113 , 8.60064926, 8.60100281, 8.60121302,
8.60133854, 8.60141277, 8.60145565, 8.60148758, 8.60151977,
8.60154455, 8.6015601 , 8.60156988, 8.60157603, 8.6015799 ,
8.60158233, 8.60158387, 8.60158483, 8.60158565])
teste alpha2_3x1_camadas
In [2]:
import mesa as ms
folder = '/media/glauffer/Data/mesa/alpha2_3x1_camadas/LOGS'
s = ms.history_data(folder)
No history.datasa file found, create new one from history.data
reading ...100%
Closing history.data tool ...
In [3]:
# extraindo os modelos, temperatura central e quantidade de zonas
model = s.get('model_number')
zones = s.get('num_zones')
ct = s.get('log_center_T')
In [12]:
import matplotlib.pyplot as plt
%matplotlib inline
plt.rcParams["figure.figsize"] = (15., 10.)
# plot do HRD
t = s.get('log_Teff')
L = s.get('log_L')
plt.plot(t, L)
plt.plot(t[154], L[154], 'r.')
plt.plot(t[2149], L[2149], 'go')
plt.gca().invert_xaxis()
plt.title(r'HRD $\alpha_{MLT}=2$', fontsize=18)
plt.xlabel(r'$\log T_{Eff}$', fontsize=16)
plt.ylabel(r'$\log L$', fontsize=16)
Out[12]:
<matplotlib.text.Text at 0x7f0b7815e190>
In [14]:
# grafico de zonas vs model
plt.plot(model, zones, 'b.')
plt.plot(model[154], zones[154], 'ro')
plt.plot(model[2149], zones[2149], 'ro')
Out[14]:
[<matplotlib.lines.Line2D at 0x7f0b788243d0>]
In [17]:
# zoom do HRD
plt.plot(t, L)
plt.plot(t[154], L[154], 'r.')
plt.plot(t[2149], L[2149], 'go')
plt.gca().invert_xaxis()
plt.title(r'HRD $\alpha_{MLT}=2$', fontsize=18)
plt.xlabel(r'$\log T_{Eff}$', fontsize=16)
plt.ylabel(r'$\log L$', fontsize=16)
plt.xlim(4, 3.5)
#plt.ylim(3.5,4.5)
plt.ylim(4,4.4)
Out[17]:
(4, 4.4)
In [ ]:
Content source: glauffer/mesa
Similar notebooks: